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0^ , 

Qs ' We study the 1/A^ expansion of a generic, strongly correlated electron model {SU{N) symmetric 

Hubbard model with U — oo and A'' degrees of freedom per lattice site) in terms of X operators. The 
leading order of the expansion describes a usual Fermi liquid with renormalized, stable particles. 
The next-to-leading order violates Luttinger's theorem if a finite convergence radius for the 1/A'' 
expansion for a fixed and non-vanishing doping away from half-filling is assumed. We find that 
the volume enclosed by the Fermi surface, is at large, but finite A'^'s and small dopings larger than 
at TV = CO. As a by-product an explicit expression for the electronic self-energy in 0{1/N) is 
given which cannot be obtained by factorization or mode-coupling assumptions but contains rather 
sophisticated vertex corrections. 
^ , PACS numbers: 74.20-z, 74.20Mn 
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I. INTRODUCTION 



Luttinger's theoremtlo is one of the most fundamental theorems in sohd state physics. It states that the volume 
V ! enclosed by the Fermi surface is independent of the strength of the interaction between electrons. One definition of 
' O ' the Fermi surface is based on the momentum distribution function n(k) of electrons. The Fermi surface is the set 
C ' of all k-points where n(k) or derivatives of it are singular. If the imaginary part of the self-energy E(k, w) of the 
one-particle Green's function G(k, lu) vanishes at lu = the Fermi surface can also be defined as the boundary of 
the set of k points satisfying ReG(k,uj = 0) > 0. If the electrons form a Fermi liquid n(k) possesses a jump and 
ImYiCkjUj) ~ uj'^ near the Fermi surface. Both definitions for the Fermi surface are then equivalent and Luttinger's 
theorem also holds. Fermi liquid behavior of the electrons, however, is not a necessary prerequisite for the validity 
of Luttinger's thcjoajm: One-dimensional metals arc Luttinger and not Fermi liquids, yet Luttinger's theorem also 
holds in this caseOu. However, there is evidence that Luttinger's theorem may be violated in strongly correlated 
^■A [ electronic systemsQQ. Based on high temperature expansions for the two-dimensional t,— J model it was concluded 
\^ . that Luttinger's theorem is violated in this case even at a doping as large as 20 per centQ. In Ref.tl finite temperature 
^^ Quantum Monte-Carlo simulations for the Hubbard model with U/t = 8 and T/t = 0.33 {U and T are the Hubbard 

0^ repulsion and the temperature, respectively) yielded a doping dependence for the volume enclosed by the Fermi surface 
^^ which was similar to that of the Hubbard I approximation and thus much larger near half-filling than that predicted 
by the Luttinger theorem. 

Theories based on perturbation expansions in njiwers of the interaction between electrons admit non-pcrturbative 
approximations which fulfill Luttinger's theorenxl. Examples of this kind are all conserving approximations. In 
i-rt strongly correlated electronic systems the interaction energy dominates over the kinetic energy and expansions around 
(~| the atomic limit. i|ij_ terms of the kinetic energy seem appropriate. Various general formulations of such expansions 

O ■ have been givenETtll. It seems, however, that so far no well-defined approximation within such an expansion has been 
O ' found iof-Pi finite dimension which explicitly obeys Luttinger's theorem. Approximate treatments based on mode- 
^ ', couplingO or factorization assumptionscJ seriously violate Luttinger's theorem. Since the employed approximations 
are uncontrolled it is not possible to judge whether these violations describe real effects or represent just artifacts. 
The use of 1/iV expansions is one way to create controlled approximations for strongly correlated systems. Formally, 



i^ • one extends the SU{2) symmetry group for the spins in the Hamiltonian to a larger symmetry group SU{N) or 
" " ' Sp{N/2) and considers 1/N as a small expansion parameter. The extension to N degrees of freedom is not unique, 
the two most popular alternatives are the slave boson or the slave fermion theories. The 1/A^ expansion can also be 
carried pjitjiicectly in terms of X operators, the analogue of the slave boson extension has been worked out in detail 
in Refs.Elltjlij. The 1/A^ expansion allows to study the validity of Luttinger's theorem as a power expansion in 1/A^. 
If this theorem holds for N — 2,3... and if the convergence radius of the 1/A^ expansion is non-zero the theorem must 
hold in each order in 1/7V. It is easy to show that the theorem holds for N ~ oo. In the following we will calculate 
the change in the Fermi surface ~ 1/N and from this the corresponding change ~ 1/N of the volume enclosed by 
the Fermi surface. In this way the validity of Luttinger's theorem for strongly correlated systems can be assessed 
under the only assumption that the 1/A^ expansion has a finite convergence radius. Explicit results are presented for 



d3/ie(13)Xe(3) + / d2dMe{123)Y{2)X^{3), (1) 



dS^ihilS)Xhi3) + / d2dMhil23)Yi2)Xhi3). (2) 



the t — J model with J = 0, the so-called t model, given, for instance, by the first term on the right-hand side of 
Eq.(l) of Ref.li3. One main ingredient of the calculation is the electron self-energy calculated in opder 1/N. Below we 
describe a more symmetric and transparent derivation of this quantity compared to that in Ref.EJ which also allows 
to determine the A'^-dependence of the various contributions in a straightforward way. We also will discuss the change 
in the n{fi) relation (n is the number of electrons per site, /i the chemical potential) if the self-energy contribution of 
0{1/N) with A^ = 2 is used. 

II. ELECTRONIC SELF-ENERGY IN 0{l/N) 

In this section we calculate the contribution of 0{l/N) to the electronic self-energy E. It represents the first 
correction beyopd the mean- field approximation for E. We consider the t — J model and use the general framework 
outlined in Ref.EiJ. Fermionic Hubbard operators, which annihilate or create an electron, will be denoted by ^e(l) and 
Xh{\-), respectively. Writing out the arguments explicitly we would have X'^'P(iiTi) and ^^^(iiri), respectively, where 
p = 1, ...^N denotes an internal index and ii,Ti a lattice site and an imaginary time, respectively. The remaining 
bosonlike Hubbard operators will often be denoted by the letter Y . Y{\) stands for ypi^i (ii^i) with pi, gi = 0, ..., TV 
with the cases pi = 0, (/i > and pi > 0, gi = excluded. It will also be often convenient to abbreviate the index 
pair JiTi by 1. All explicitly given formulas in the following refer to the t — J model with J = 0. Using the rules 
discussed in RefJl3 it is straightforward to generate from them the corresponding contributions due to a spin-spin or 
a charge-charge interaction term. 

We start from the equations of motion for X operators, 

dX,{l) 
dn 

dX^jl) 

The functions /i and t in Eqs. (M and (0) are defined by 

Me(12)=/i<5(I-2)5,,,„ (3) 

^lh{l2)^-^lS{l-2)Sp,p„ (4) 

ie(123) = — (5(1 - 2)[(5p2 0'5g2 0^gig3 + ^<!iq2^P2q3]^ (5) 

t/i(123) = — (5(1 - 2)[6p^oSq2o5pip3 + Sp.pjq^ps]- (6) 

The electronic Green's function G is defined by 

G(12) = -(rXe(l)X„(2)), (7) 

where T is the time ordering operator. Dyson's equation for G reads: 

G-i(12)= I d3Q~\l3) [Gq\32) + ne{32) - J:e{32)] . (8) 

The free Green's function Gq is given by 

d 

Go ^(12) = -^(1 - 2)[Sp^oSp^oSqiq2 +Sqi05q20SpiP2]-^- (9) 

There is no need to specify any electron or hole index in Gq^ since Gq ^ is identical in both cases. Q is the spectral 
weight given by the expectation value of the equal-time anticommutator of the two fermionic operators of G, i.e., 
writing out explicitly the internal indices, by 



1 2 



Qi T o ) = '^(i - my""a))Sp, + (Y'^^'im- m 



The self-energy Se(12) is given by 

Ee(12) = / d3d4d5ie(134)(Ty(3)X(4)X(5))G-i(52). (11) 



Using Eq. (Q) we apply the operator G ^ in Eq. (|ri| ) to the left and act with it on the time-ordered product of Hubbard 
operators. Using Eqs.(E])-(0) and evaluating equal-time anticommutators we obtain 

I]e(12) = - /d3d4d5t(1345)(y(3)y(4))Q"i(52)-|- /d3d4d5A(1345)(X,,(3)Xe(4)>g-^(52) 

d3d4:d5d6d7te{l3i)th{756){TY{3)Y{5)Xe{4:)Xh{6))Q-\72), (12) 

with 

t(1342) = — S(l - 3)S{2 - i)[Sq^pJpsoSq^aSp^aSq^o + Sq^qJp^pJp^oSq^a 

'^^q-^q^Op.^p^Op^QOq^O + "gigs "P3g4''p4P2 J J \^'^) 



and 



N 



A(1342) = ^^^<5(I - 2)^(1 - m,^pJps,. - K.Jp^P.]- (14) 



In Eq.(|12D we used the fact that reducible self-energy contributions cancel each other and thus can be dropped. Since 
this property of self-energy diagrams is not completely obvious in our case we give a proof of it in appendix A. The 
first two terms in Eq. ([l2|) are frequency- independent and thus contribute only to the dispersion of quasi-particles. 
The third term in Eq.(|2p depends on the frequency and determines the damping of the quasi-particles. 

It is convenient to introduce an external source field K coupled to bosonic Hubbard operators for the evaluation of 
the expectation values in Eq.(|lJ) similar as in Ref.EJ. The self-energy assumes then the form 

Se(12) = - /d3te(132)(r(3))- f d3dM5t{13A5)^^-^Q-\52) 
J J dK{3) 

d3d4d5A(1345)(X(3)X(4))g-i(52) 

d3d4d5d6d7d8d9ie(134)t,i(956)G(47)^^^^^^G(86)Q"V92) 

oK(6) 

d3d4d5d6d7d8d9dl0dlUe(134)tft(1156)G(47)r(78;5)G(89)r(910;3)G(106)Q-^(112) 

d3d4d5d6d7te(134)t,,(756)^^^^G(46)Q"\72). (15) 

oil (3) 

r is a vertex function defined by 

We tacitly assumed in Eq.(|l5[), and will do so also in the following, that the external source K is put to zero once the 
functional derivatives have been carried out. The evaluation of the fourth term in Eq.(nq) is somewhat more involved 
than that of the other terms and is therefore deferred to appendix B. In the following we often will abbreviate this 
term by ST.. Summing over the internal indices or exhibiting them explicitly Eq.(n5^ becomes 

S-(I2) = -fc^ [{Y°°) + {Vn] +iN- l)fc^5(3i+)<5(i - 2) 

. ^(i-3) ,, t(4-2) t(I-2) 

^ TV ^^^ N N[{YOO) + {Y'^'^)]' 



j^yaio-1^2)) 6{Y'"'{2)) (5(y'"i'^i(2)) 5{Y'"'^{2)) 



t(I-3) 



5K°°{\) 



SK'"'{1) 



N 



5(34)5(56)5(78) 



t{8 - 2) 

TV 



[(r0") + (y^^}]2[r^-(45;2)r; 



00 



(67; I) 



TZl (45; 2)r^^^^ (67; I) - r^^(45; 2)r--(67; T) 



+r^^^U45;2)r^J-(67;I) 

+r^r (45; 2)CJ^(67; I) + T^^J (45; 2)r^r (67; I)] + <5I]^^(i2 



(17) 



In Eq.(0) we have renamed the internal indices p,pi... by a,(Ji... where u, ui... assume only values between 1 and N. 
We also used the fact that in equilibirum E and g are diagonal and independent of spin indices, Y is diagonal and in 
r either the rows or the columns have the same spin indices (the upper two indices in F result from the two fermionic, 
the lower two indices from the bosonic Hubbard operators) . Spin indices on the right-hand side of Eq. ( [l7|) which do 
not appear on the left-hand side of this equation are summed over. In accordance with that we explicitly show the 
index a, for instance, at the self-energy, in spite of the fact that S is independent of spin labels in equilibrium. In 
Eq.(pT[) we also have introduced a normalized Green's function by means of 



G(12)= /d35e(13)g(32)= /d3Q(13)5/,(32). 



(18) 



In equilibirum we have ge — gii so we simply can use the letter g in Eq.(|l7|) for convenience. Let us introduce 
normalized vertex functions by 



7e(12;3) 



'^.9e'Hl2) 

5K{i) ' 



(19) 



7/.(i2;3) = 



^.g,"^(12) 
6K{-i) ' 



(20) 



Unlike ge and 5/1, 7e is not equal to 7/1 in equilibrium. The previously defined vertex function T can now be expressed 
as 

r(12;3)= /d4Q-i(14)7e(42;3)-H j dA%^g-\A2) 



5K{3) 



d47^(14;3)Q-i(42)-f / ^45^^(14) 



/g'^(42) 
SK{3) ■ 



(21) 



We stress that the two expressions for F in terms of e- or h- quantities are identical, so that we may take the more 
convenient one of the two whenever we use Eq.(|2l]). 

Let us introduce some short notations that will be useful later: 



n 



(TiO"2 0"3(T4 






(22) 



(7 



CTi £72 (73 (74 



2- 3 ,=KIZ{12;5), 



(23) 



Hi ^ 

The 1/N expansion can be simplified by noting that 



2 . 0'3'74 

2 ' 3 






K=0/ 



7^3^^! (12; 3). 



<0{A 



a/5 I 



\K=Ol ' 



(24) 



(25) 



where A°'^ is a functional of 5. 0(...) on the left- and right-hand sides of Eq.(25) denotes always the largest order 
occurring in one of the tensor elements. The functional derivative on the left-hand side of Eq.(|25|) acts on a Green's 
function g'^^'^'^ and can be written as 



omitting space and time labels, zeros in fermionic index pairs and taking finally all quantities in equilibrium. If 
f3 = a4 in K the number of internal sums on the two sides of Eq.( |25|) is equal. Since the orders 0{g) and 0(7) are 
at most 1 Eq.(p5|) holds. If 0-3 ^ a4 in K the left-hand side of Eq.(p5|) may contain first one sum more than the 
right-hand side of this equation. However, because of Eq.(|26|) and because of the selection rules for 7, this additional 
sum does not survive when finally going to equilibrium so that Eq.(^5|) again holds. Using Eq.(p5|) the expansion of 
various quantities in powers of 1/iV read 

E(i2)=^E.(i2), (27) 

i=0 

00 

g(l2) = ^5,(12), (28) 

j=0 



(roo) = ^ (rO")„ (29) 



(y-}=^(r--)„ (30) 

i=0 



n(2;T)=£n,(2;I), (31) 



r(l2;3)=^r,(T2;3), (32) 



7(T2;3) = ^7«a2;3). (33) 

1=0 



Carrying out the 1/A^ expansion, we find for the leading order of the self-energy: 



/yOO] 
Eq.(^4[) reads in Fourier space 



S^-(i2) = - ^^ ^ ' t{\ - 2) + t{\ - 3).9o(3l+)<5(I - 2). (34) 



I]o(k) = -^^^i(k) + Ao, (35) 

where 

Ao-^t(k)go(fc)e^""°^ (36) 



and 



lujn + ((y"'J)-i/iV)i(k) -f ^ - Ao 



In Eq.(36) and in the following equations the short notation ^j, stands for {T /V) X]k ^j where T is the temperature 



and V the volume of the crystal. 



For the order 0{1/N) of S we obtain: 



£^'"(12) 



i(i-2)r,^oo 



N 



[(r"")o + (r--)o] + t(I-3) 



5i(31+) 



+^ 3o(31+) 



N 



5{l - 2) 



,i(l-2) i(I-3) 



+ ( /.Vnnx - ' .. go 

i(l-3) 



(34M4 - 2)) [no [-,'1) +no ("2^; Y) - no (Y; ^ 



' 5o(34)5o(56)go(78)i(8-2)(y"0)2_, x 



TV 



y-, ( a. a 00\j-i ( a. a aa\ j^ f a a era \ y-, fa a 00 

.^H45'2J^H67' lJ+^H45' 2J^H67' T 

T^ /aa ctctAt^ / a a aa\ j^ ( a a aa\ j^ ( a a <t<t\1 . cv^(T(T /Ton 

+^H45' 2J^H67' lJ+^H45' 2J^H67' iJJ+'^^i ^^^^^ 
We are now evaluating analytically all the unknown quantities in Eq.(pq) using the relatiouEiJ 



(>^^^(i))iK=o = .Hii+J-nii'i 



(r"0} 



1 1- 






and the constraint 



N 



yOO _|_ y^ yaa 



N 



a=\ 



Performing explicitely the 1/7V expansion, we obtain 



N 



{Y^'^')^, = --N{Yno, 



(y"")o = -N{Y^^)i, 



{rno = 90 (n+) , 



{¥--), =g, (11+) + 



(roo)_i 



-no (Y;°^) +.9^(11+) -5o(n+) 



no("2";"j)=-5o(23)7e,o(g|;™)5o(42+), 

no ("2"; Y) - -.9o(23)7e,o (ll; Y) .90(42+), 

no (Y ; 7) = -5o(23)7e,o (ll; Y) 5o(42+), 

7e,o ( ^ ^ ; I) = -<5(1 - 2)5(1 - 3) + [J(I - 2)t(T - 4) + i(I - 2),5(I - 4)] 
xffo(45)7e,o(gg;°|)5o(6l+), 

7,^0 ( J ^; 1) - -<5(I - mi - 3) +go(24)7M (^ gi "j) 5o(56+) 
X [,5(1 - 2)i(6 - 2) + t(T - 2)(5(6 - 2)], 



7e,0 



( J ^; 7) ='^(1-2)5(1-3), 



(38) 
(39) 
(40) 

(41) 

(42) 
(43) 

(44) 
(45) 
(46) 
(47) 

(48) 

(49) 
(50) 



^ (aa 00\ 1 Zero- 00 \ , ^ ^ ( aa 00 \ -1/THn 

1 (a a 00\ , -^ -1/ToMT ( '^'^ 

- (yooyi^^M ii 2' 3 J + (^^ooyr^so (i2)no y ^ ; 



(51) 



00 

37' 



HI 2' 3 J " (yoo)_J'=^Hl 

( a. OL aa \ -^ * 

VT2' 3/ 



Cr OL _ (7Ct 

2' 3 



(r"0)2_^' 



olol aa \ — 1 / 



ocn^no(7;7)<7o-^(12) = 0. 



Furthermore 51 is related to go by 



gi(l2)=.9o(l3)Si(34)go(42), 



(52) 
(53) 
(54) 

(55) 



which yields a self-consistent equation for gi, or, equivalently, for Si. 
Using Eqs.(|9[]5|) we obtain for Y^l" : 



S?'"(12) 



[(y°°)o + (y-)o] 

N 



i(l-2) + (5(l-2)t(l-3) 



ffi(3r 



.9o(3i+) 



N 



i(l-2) i(I-3) 
_(FOO)_i ]V 

i(I - 3) 



5o(34)t(4-2) 



no(-;-)-n„(-;Y 



^ _ 
_N _ 

^(i-3) 

TV 



-50 



(34)7e,o(5g;|)ffo(5l)ffo(i6)i(6-2) 
5o(32).9o(24)7.,o (| ^ ; °°) 5o(56)t(6 - 2) 
5o(34)/i(45; 12)5o(56)5o(i2)ffo(27)i(6 - 2) 



(56) 



In Eq. ( pq ) we have used the explicit expression for the term bTj^" (Vl) calculated in appendix B where the function 
/i(45; 12) is also defined. 

The expectation value of the particle number operator with spin a \s n„ — (y^'^)o + ^'"')\ — nga- + i^icn where 
the_expectation va_lu_es are_given by Eqs.(||) and (||). Writing 7(l2; 3) = 7(1 - 2, 1 - 3), n(l2) = n(I - 2) and 
/i(12;34) = h{A — 2,4 — 1,1 — 3), performing Fourier transforms, and adding the contributions of 0(1) and 0{1/N) 
the self-energy becomes 



A^- 1 



N 



N ~l 
-1 H -j;r—2no^ + 2niCT 



2 



Nil - 2noa) 



Y, t{k + q)An{q) - ^ E *'(k + q)5o(fc + q)AU{q) 



--^ X! ^(^ ~ ^i^^^P + ^)9o{k - q)go(p + q)9o{p)['ie^o{p, q) + ihfi{p + q^ -q)] 



P,9 



-j^^iiv^goik + q))gQ(p + q)ga{p)lefl{P-.q)^ 



c{q) 



b{q)' 



where AU{q) = U^giq) - n-^(g) and 

ai =^i(k)gi(fc)e-"°" =5]t(k).92(fc)Ei(fc). 



(57) 



(58) 



Here we wrote the second index pair of Il{q) as an subscript to simplify the notation. Explicit expressions for the 
functions 7, 11, ft, are 



l + b{q)-aiq)t{k) ^^^^ 

(60) 

n^S(g)=a(q), (62) 



Ie,0\'^,H) — 


[l + b{q)][l + b{-q)]-a{q)c{q) 


-fh,o{k,q) = 


l + 6(-9)-a(g)t(k + q) 


[l + b{q)][l + h{-q)]-a{q)c{q) 


r\'^'^(r,\ — - 


a{q) 



h{ki,q,k2) =7e^o(fci, 9)90(^2) 
and for the susceptibilities a, b, alii 






I + K9) 



(63) 



a{q) = -Y,9o{k)go{k + q), (64) 

k 

b{q) = -Y.mgoik)go{k + q), (65) 

fc 

c{q) = - 51 *(k)i(k + q)5o(fc)5o(fc + ?)• (66) 

fc 

Note that je.oip, q) = 7/!.,o(p + <?, -<?)• 

Finally, the relation between the particle number n^ and the chemical potential becomes, using Eq.(E3), 

"' - "- (' - ^T^) - W^) ? "" "> + *"■'■ '"' 

where 5nia is defined by 

5n,, = ^ gi(fc)e*"""^ = ^ 5o'(fc)Sr(fc)- (68) 

fc fc 

Eq.(p7|) gives an explicit expression for the self-energy taking into account the orders 0(1) and 0{1/N) of the 1/A^ 
expansion. The first two terms represent renormalizations of the chemical potential. The third and fourth terms 
are frequency-independent and thus contribute only to the dispersion of quasi-particles. Terms 5-7 are frequency- 
dependent, contribute to the damping of quasi-particles, and cannot simply be obtained by a factorization or mode- 
coupling assumption applied to Eq.(0). 

Eq .(|67|) represents the relation between the particle number and the chemical potential. As can be seen from 
Eqs.([43D and ( p4| ) n^ is not only determined by the zeroth and first orders of the normalized Green's function g but by 
additional terms. The appearance of such additional terms is not unexpected because n„ is given by the normalized 
Green's function g(ll+) at TV = 00 and by the Green's function G'(ll+) at TV = 2 so that the \/N expansion must 
interpolate between these two cases. 

III. VIOLATION OF LUTTINGER'S THEOREM IN 0(1/N) 

Eq. (p^ shows that the imaginary part of the self-energy vanishes at the chemical potential. There is thus a Fermi 
surface {ki?} determined by the condition 

M = S(kf,c^ = 0,/i), (69) 



where the chemical potential /i is determined by the particle number n = Nria per site, ria does not depend on 
the spin direction a in the normal state, of course. Nevertheless, it is convenient to denote the particle number per 



spin direction by n^ to distinguish it from the total particle number n. In Eq.(69) we have written explicitly the 
dependence of S on the wave vector k, the frequency co and the chemical potential /i and omitted unnecessary spin 
labels. 

In the unconstrained, free case the self-energy S in Eq.(^9|) should be replaced by —t(k)/N. Each state k can be 
occupied by N particles. The set of Hamiltonians with different N's used for the 1/A^ expansion should correspond to 
the same Ua-. As a result the free Fermi surface {k^?} is independent of N though the corresponding chemical potential 
fj, and the electron dispersion depend on N. Turning on the interaction between electrons the Fermi surface depends 
in general on TV. Luttinger's theorem states, however, that the volume enclosed by the Fermi surface is independent 
of N and also of the interaction strength. 

In the constrained case one expands for a fixed Ua- the quantities /i, S, and hp in powers of 1/A^ similar as in 
Eqs.@-(|3|). Using Eqs.Q and (|l]) the zeroth order of Eq.® is 



Eqs.(|3) and (E|) yield at T = 



Mo = -(2 ^"-ocr)i(ko_F) + Ao(mo)- (70) 



noa=X!^(^o-Ao + (--no^)i(k)). (71) 



For a fixed noa Eq.(O) determines /ig — Aq and Eq.([70|) ko_F. These equations are the same as in the unconstrained, 
free case if we choose (i — noCT)t(k) for the dispersion and renormalize the chemical potential appropriately. This 
means that at A'^ = 00 the Fermi surface is identical to that of the free, unconstrained case and that, in particular, 
Luttinger's theorem holds. This result is not very surprising because at iV = co the particles are renormalized but do 
not interact with each other. 

The first-order of Eq.(p9) in 1/A^ is given by 

aEo(koF,/xo) , 9So(koF,Mo) , , ^ ,, r, ^ ,„„^ 

A*i = 5 fJ-i H ^1 kiF + Si(koF, w = 0, fio). (72) 

Ofio dkoF 

Here kiF is assumed to be parallel to the vector d'So/dhop with kpF as origin. Because Sq is independent of frequency 
we have dropped the frequency argument in Eq for convenience. The change /ii in the chemical potential is to be 
calculated for a fixed n„, i.e., from 



''(7^ 



{dnoa/d^io) ■ Ml + niaino) = 0. (73) 

Eq.([72[) must hold for any point kpF on the Fermi surface. This equation thus determines kiF as a function of kpF 
and thus the change of the Fermi surface in 0{l/N). Note that due to the last term in Eq.(|7^) the shape of the Fermi 
surface will in general change by the interaction between the particles. 

Eq.([72|) can be simplified in the following way. S depends on fi only via /x — A where A is again a function of /x. 
Defining a renormalized chemical potential /i by 

A = M - A, (74) 



we can drop A in S everywhere replacing /x by jj. Eq.([73|) also holds if /xg and /xi are replaced by /xq and /xi, respectively, 
yielding 

(dnoa/dfio) ■ Ml + nia{flo) = 0. (75) 

It is then easy to see that the first term in Eq.([72[), written in terms of /i, cancels the term nicri(k) in Eq.(p7[). Eq.([7^) 
thus becomes 

9So(koF,Mo) , , fi n r, ~ \ ('7a\ 

Ml = TT, kiF + 2]i(koF,w = 0,Mo), (76) 

Ckof 

where Si is given by Eq.(^) with all terms in the first line on the right-hand side dropped except for the term 
— nocrt(k)/A^. According to Eq.(fl]) the derivative dnoa/dfio gives rise to two contributions, one in which the derivative 
acts on /xq and one, where it acts on 71o<t- We also split riia, given by Eq. (p7\) , into three contributions. 



ni„ ^ ni„ + Snlcr + 5ni„ , (77) 

with 



ni 



' N{1 - 2noa) 



(l + Y^U^oSil)), (78) 

Sh^, = }_^gm^rik), (79) 



k 



5ni, = Y,gl{k)ni,t{\^). (80) 



Using these results Eq.(|73|) can be rewritten as 



„- i^u: ■tJ:i+y^Ql{k)( ^r^)t{k)fii+hi„ + Siller + Snicr = 0, (81) 

dfj.0 ^ V d^ia / 

wherejex means that the derivative should be taken only with respect to the explicit dependence on JIq. ^From Eq.(73) 
and ( pO| ) follows that the second and fifth terms on the left-hand side of Eq.(8^) cancel each other. Solving Eq.(80|) 
for fii we insert its solution into Eq.(|76|) and average the resulting equation over the zeroth-order Fermi surface and 
obtain 



fiia - 5ni„ = ^(5(ek ~ Ao)-^ • ki_F +^(5(ek -/io)Si(k,w = 0,/2o), (82) 

k k 



where et denote the zeroth-order one-particle energies. The first term on the right-hand side of Eq. (p2f) is equal to the 
Fermi surface integral / dS-p ■ kii? — V^ , where Sp is a directed element of the Fermi surface and Vf is the change 
in 0{\/N) of the volume enclosed by the Fermi surface. Eq.(p2|) becomes finally 

V^ = -fii, -Y,gl{k){ti{k,ilo) - Si(k,c^ = 0,Mo)). (83) 

k 

The first term on the right-hand side of Eq. ( p3|) is due to the unusual relation between the particle occupation 
ricr and the chemical potential jl specified in Eq. (|67[) . In Fermi liquid theory ria is given by the one-particle Green's 
function integrated over all momenta and over frequency up to ft. Such a relation holds in our case only at iV = oo 
and N = 2 with g and G to be used as Green's function, respecively The 1/A^ expansion must interpolate between 
these two cases which is the origin of the more involved nature of Eq.(|67|). Thus the first term in Eq.(^3|) is intimately 
linked to the 1/7V expansion and of fundamental nature. In Fermi liquid theory it is, of course, absent. 

In the second term on the right-hand side of Eq.(|83|) any frequency-independent contribution to Si drops out. 
This term has the same form as the term in Fermi liquid theory where the product of the frequency derivative of Jie 
self-energy and the Green's function is integrated over (second term in the square bracket in Eq.(19.12) in Ref.EZi). 
The only difference is that in our case the involved Green's function is the zeroth-order Green's function g^ and E 
is a functional of gQ. This term vanishes in Fermi liquids for general reasons yielding Luttinger's theorem. Though 
we are presently unable to prove something similar for the second term in Eq.(|83|) it is easy to see that it cannot 
cancel in general the first term on the right-hand side of Eq.(B3). ^From Eq.(|7q) follows that ni diverges in the 
limit noa — > 1/2. In contrast to that the second term in Eq.(S3) is clearly regular in that limit. /^From this one 
concludes that there exists at least a finite, A^-independent interval in doping near half-filling where the first term in 
Eq. (p3[) dominates and makes V^f non-zero. Assume now that V^ is zero in an arbitrarily small interval of the doping 
6 — 1 — 2norj. This implies that the regular part of V^ is equal to A/ 6 with some constant A throughout this interval. 
Performing an analytic continuation in the variable 6 in Ei and Vf^ shows that the regular part should be given by 
A/ 6 in the whole interval < S < 1 ii there is no phase transition present making the interchange of sums over k and 
the continuation in S impossible. The regular part i n V^ would then develop a singularity at (5 = in contradiction 
to our assumptions. Since detailed calculations showlla^tj that at [/ = oo and large A^'s there are no phase transitions 
one concludes that the Luttinger theorem must be violated for dopings throughout the entire interval < S < 1. 

Fig. 1 shows the dependence of the sum of no and ni as a function of no for A^ = oo (dashed line) and N = 16, 8, 4, 
and 2 (solid lines), ni is always negative which means that for a given uq or /xq less particles can be accomodated at 
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FIG. 1. Sum of the 0(1) and 0{1/N) contributions no and n\, respectively, to the site occupancy as a function of no for 
different N^s. Only terms which diverge near half-filling haven been considered in n^. 

a site if electronic correlations are taken into account compared to the uncorrelated case. This effect is very small at 
small no's but increases monotonously towards half- filling where ni diverges. For A^ = 2 a maximal site occupation 
of only about 0.4 is possible. For larger /io the n(/io) relation has a maximum and then a negative slope signalizing 
an unstable homogenous state. This result clearly shows that a calculation which includes only 0(1) and 0{1/N) 
contributions cannot be used to describe the physical case of A^ = 2 near half-filling. With increasing TV the maximal 
site occupancy increases monotonously towards 1 and the maximum of the curves is shifted towards ng = 1. 

Fig. 2 shows the total volume V^ enclosed by the Fermi surface and calculated in 0{1) and 0{1/N) keeping the 
singular first term in Eg. (p3[) . The dashed line corresponds to A^ = oo and represents exactly the volume predicted 
by the Luttinger theorem. The solid lines describe the cases N = 16,8,4, and 2. The Figure shows that electronic 
correlations, taken into account in 0{1/N), always increase V^ . This increase is small at small electron concentrations 
but becomes larger and larger at small dopings. Our results are asymptotically exact at sufficiently large iV's and 
in a finite, iV-independent doping interval near half-filling. As argued above this implies that Luttinger's theorem is 
violated throughout the entire doping interval < S < 1. This means that the solid curves for large A^'s are never 
coincident with the dashed line in Fig. 2 and that these lines are at least qualitatively a good approximation for all 
dopings. Fig. 2 illustrates also the fact that increasing electronic correlations increases the Luttinger volume V^, 
i.e., V^ is positive. The solid lines in Fig. 2 diverge towards half-filling due to the 1/A^ expansion of the reciprocal 
spectral weight Q~^. The convergence radius of this expansion is determined by 1/A^ < S/hq and thus becomes very 
small at small dopings. This means that for a fixed A^ the power expansion in 1/N breaks down at small dopings. 
The larger A'^ is the smaller is the critical value for the doping where this expansion breaks down. 




1-no 

FIG. 2. Sum of the 0(1) and 0{1/N) contributions to the Luttinger volume V^ as a function of the doping 5 — 1 ~ no 
keeping only divergent terms in the 0{1/N) contribution. The dashed line corresponds to A'^ = cx) and obeys Luttinger's 
theorem. 
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IV. CONCLUSION 

Using a 1/N expansion we have derived explicit expressions for the 0(1) and 0{1/N) contributions, So and Si, 
respectively, for the self-energy of a SU(N) symmetric Hubbard model with N degrees of freedom and U = oo. Using 
previously derived rules these results can immediately be generalized to the case of a finite, but large U . We find 
that the frequency-dependent part of S is regular in the doping i5, whereas the frequency-independent part develops 
a singularity in 5 due to the expansion of the inverse spectral weight Q~^ in powers of 1/A'^. This means that near 
half-filling large A'^'s are required to obtain quantitative results. As an example we find that the n(/i) curve, calculated 
with So and Si and iV = 2, becomes unreliable already near and below a doping 5 ^ 0.20. Our derivation of S 
also shows that Si involves rather sophisticated vertex corrections and that usual many-body assumptions (neglect 
of vertex corrections, lowest-order iterations of the functional equations for the self-energy and the vertex) cannot be 
used in order not to miss cancellation effects. These features are connected to the fact that the hopping t cannot be 
considered as a small parameter in the perturbation expansion in the kinetic energy of the electrons, at least in the 
metallic state. This is related to the fact that t appears not only in the numerator but also in the denominator of the 
perturbation series because the dispersion of the one-particle energies is also determined by t. 

Using the above results we have studied the 1/A^ expansion for the Luttinger volume V^ — Vq + Vf + ... where 
V^ denotes the volume enclosed by the Fermi surface. The 0{\) contribution V^[ obeys Luttinger's theorem but this 
theorem is violated already in 0{\/N). We have shown that V^ ^ throughout the entire doping interval < 5 < 1, 
if the normal state is stable, and that V^ > at least in a finite, A^-independent interval near half-filling. Keeping 
only terms which are leading at small dopings our calculations suggest that V^ > for all dopings and that V^ 
becomes very small at large dopings. There are also good reasons to believe that the normal state to which our 
discussions is restricted, can be considered as the stable state at large A^'s throughout the intesyal < S < 1. The 
largest transition temperatures to superconducting states behave roughly as ~ C/~^ at large C/'st3. Strictly sneaking, 
this is not entirely true because extremely weak superconducting instabilities have been found even for U = cxxj. The 
corresponding transition temperatures are, however, astronomically small so that these instabilities can be neglected 
for our purposes. In a similar way all the other instabilitiesllj of the t — J model at large A^'s such as flux phases, 
bond-order or charge density waves, vanish in the limit J — > or f/ — > co. This means that we have indeed dealt with 
the most stable ground state of our system. 

Finally we want to poJJit out that our results agree with two recent numerical studies of the Luttinger volume in the 
Hubbard model. In Reffl Quantum Monte-Carlo simulations for the 2D Hubbard model with U/t — 8 and T/t — 0.33 
have been preformed and an increase in Vf beyond the value predicted by the Luttinger theorem has been found for 
S < 0.2. Since at these high temperatures structures due to the quasi-particle peak and the incoherent background 
cannot be distinguished the interpretation of the results is not unique: a) Vp may be the volume enclosed by the true 
Fermi surface or b) Vp is a volume enclosing a "spectral weight Fermi surface" and thus be related to the spectral 
distribution of the .whole Green's function. Similar results, subject to a similar uncertainty in the interpretation, have 
been found in Ref.u using high-temperature expansions and considering the momentum distribution function. In our 
case it is clear that only the case a) can apply and that Vp denotes the volume enclosed by a constant energy contour 
of quasi-particles. 
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APPENDIX A: CANCELLATION OF REDUCIBLE SELF-ENERGY CONTRIBUTIONS 



In Eq. (15) we have omitted the foUowmg reducible self-energy contributions: 

E^^''(12) = - f dl'd2'd3d3'dAte{132'){Y{3))G{2'l'){Yi3'))t,,{-i3'l')Q-\A2) 



dl'd2'd3d3'd4'd5'fi6ie(132')(r(3)}G(2'l')r(l'4';3')G'(4'5')i/i(63'5')Q"^(62) 
+ f dl'd2'd3d3'dM5d6te{13i)G{A5)r{52';3)G{2'l'){Y{3'))th{63'l')Q-\62) 
dl'd2'd3d3'd4d4'd5d5'd6ie(134)G(45)r(52';3)G(2'l')r(l'4';3')G'(4'5')i/»(63'5')(3"^(62) 
dl'd2'd3d4(i5d6ie(134)G(45)r(52';3)G(2'6)Q"i(6l')Se(l'2) 



dl'd2'd3d4ie(132')(F(3))G(2'4)Q-i(4l')Se(l'2). (Al) 

In equilibrium the factors Q^^{12) yield a delta function (5(1 — 2) and a constant factor which can be taken outside of 
the integrals. Writing the expectatin value of three Hubbard operators in Eq. (|l l[) in terms of a functional derivative 
and using Eq.(p^) we obtain from Eq.dnl) 



Ee(12) = - / d3te(132)(r(3)) + / d3d4(i5ie(134)G(45)r(52; 3). (A2) 

An alternative procedure to calculate G uses the equation of motion Eq.(0) instead of Eq.([l|). Eqs.(p|) is then replaced 
by the equivalent equation 

G-i(12) - / d3[Go\l3) - A^„(13) + E„(13)]g-i(32), (A3) 



and Eq.(|A2|) by 

T.h{12) = - I d3th{231){Y{3)) + I d3d4d5r(14; 3)G(45)tft(235). (A4) 



In view of Eqs.(|) and (|) Eqs.(|) and (^) yield E;,(12) = -Ee(12). Using this identity as weU as Eqs.(^) and 
(A4) Eq.(|Al[) can be written as 

E^'='^(12) = -g-i f dl'd2'I]e(12')G(2'l')S/,(l'2) - Q-^ f dl'd2'I]e(12')G(2'l')Se(l'2) = 0. (A5) 

We thus have shown that the reducible self-energy terms cancel each other and do not contribute to the self-energy. 
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APPENDIX B: EVALUATION OF ST, 



Let us define a function F by 



F(12,34) 



d^G-^12) _ (5r(12;3) _ (5r(12;4) 

SK{3) ■ 



6K{3)SK{A) SK{A) 
The equilibirum contribution ST. can then be written as 



(Bl) 



^£'"'"(12) 



d3dM5d6^-^^-^g{M)g{m)*-^:^—^ x [{Y"°) + {Y"")] 



F 



a a- 00 00 



45 2 



CTi crcri 00 1 i^ 

45' 2 i '+^ 



i)-^a 



/(Ti a ^ 00 (Tio-\ f/'^i °'' 

V4 5'2i/ V4i 



2 CJ"CJ"2 cri(T 

5 ' 2 1 



(B2) 



where we have written out the spin labels and summation over indices which do not appear on the left-hand side of 
this equation is assumed. It is also convenient to introduce the function 



/(12,34)== 



SK{3)SK{4:)' 



(B3) 



F can then be expressed as 



i^(12; 34) = Q-\l)f{U; 34) + ^-%^^,{U; 4) 



SQ-\l] 
SK{4) 



7e(12;3) 



SK{3) 
SK{3)6K{4) 



.9e-'(12). 



(B4) 



Next we expand the functions F, f in powers of 1/N similar as in Eqs.(p^)-(p3). Considering the 0{1/N) contri- 
bution of ST,, i.e., ST,i, the first term in the square bracket in Eq.(B2) drops out and we obtain 



ST'mu) 



d3d4d5d6^^i— ^go(34)5o(56)<(6 - 2){Y°")_^ x 

j-i / a a era ^^\ \_ rri f ^ ^ 00 a(T\ ^ r f ai a aii 

H45' 2 I J H45' 2 1 J (ro")_/U 5' i 



1(7 era 

2 



(B5) 



In the last term in the square bracket in Eq.(B5) the 0(1) contribution is to be taken and there is a sum over ai 
whereas a is a fixed generic index a ^ a. Straightforward considerations show that 



r f a a era 00 A r / a cr 00 aa \ ^ 

■^Hl2' 3 4J^^Ht2'3 a)^^' 



(B6) 



-Fi 



{(TO. 

V12 



aa 00 

2' 3 4 



^1 



CT Q _ 00 aa. 

i 2' 4 3 



^^no(^^;°|)<5(I-2)<5(I-3). 



(B7) 



For example, consider Eq. (|Bq ). We rewrite the functional derivative 5/5K{ ^ ) as a derivative with respect to g 
and obtain in leading order 



^°\ir 3 



SK{°i) 



J2 / d6d7 



'^7(11; 



6 7/ 



■ 50 (64)50 (47), 



(B8) 



where all quantities are to be taken in equilibrium once the derivatives have been carried out. The 0(1) term of 7 
is independent of g and thus drops out. Since in the generic case f3 ^ a, f3 ^ a, one internal sum in the 0{\/N) 
contribution for 7 is removed by the derivative SjSg (f 1) but the sum over (3 gives a factor TV. This means that the 
right-hand side of Eq.(p^) is of order 0{1/N) and yields Eq.(p^). 

Finally, the last term in the square bracket of Eq. JBq ) can be evaluated as follows. First we rewrite it in a similar 
way as /o in Eq.(B8) and obtain 



•^ V 4 5' i 2 



5-1 \ 



4 5' 1 



6K\ 



d6d7 



(57 I 



4 5' 1 



.50(62)50(27). 



(B9) 
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The quantity S^/Sg in Eq.(B9) is very similar to the quantity g'^^^ in Eq.(34) of Ref]l3 and can also be calculated in 
a similar way. One obtains 



Sj 



4 5' 1 



S{6-l)-h{A5;17), 



where h satisfies the equation 

/i(45; 17) = / d3i(4 - 3)7e,o(75; 4)5o(3 - T) 

d3d6t{4: ~ 3)go(3 - 6)go(4 - l)/i(65;47). 



Writing /i(45; 17) — h{7 — 5, 7 — 4, 4 — 1), and performing Fourier transformations, Eq. ( Bll ) becomes 
h{ki,q, k2) = t(k2)5o(fc2)7e,o(fci, g) + 30(^2) ^ t{p)go{p)h{ki,q,p + q), 



(BIO) 



(Bll) 



(B12) 



with the solution 



h{ki,q,k2) = 7e, 0(^1,9)50(^2) 



i(k2) 



c{q) 



1 + b{q) 



Then, inserting Eqs.(B7) and (P38D into (B5), and making use of Eqs.(|B9|)-( Bll ), we obtain 

2 



6^rik) - -^TSjyT^^o J2 *(k + q).9o(fc + 9)n^o (g) 



-j^'^t{p)9o{k + q})go{p + q)9o{phefi{p, q)^ 



ciq) 



Hqy 



(B13) 



(B14) 



The first term in Eq. (B14) cancels with a similar contribution in the total self-energy (38) coming from the 
terms Fi (| | ; "l) Fi ( g ^; y) and Ti (| | ; °'2° ) Ti ( g 7 ; '\") once the corresponding explicit expressions Eq.(|5|) for 
the charge vertex are used. 
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FIG. 3. Sum of the 0(1) and 0{1/N) contributions no and n\, respectively, to the site occupancy as a function of no for 
different A'''s. Only terms which diverge near half-filling haven been considered in n\. 



17 




FIG. 4. Sum of the 0(1) and 0{l/N) contributions to tiie Luttinger volume V as a function of the doping (5 = 1 



no 



keeping only divergent terms in the 0{\/N) contribution. The dashed line corresponds to A'^ 
theorem. 



oo and obeys Luttinger's 
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